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Abstract 

We study the non-nearest-neighbor interaction effect in 1-D spin-1/2 chain model. In many 
previous schemes this long-range coupling is omitted because of its relative weak strength compared 
with the nearest-neighbor coupling. We show that the quantum gate deviation induced by the 
omitted long-range interaction depends on not only its strength but also the scale of the system. 
This implies that omitting the long-range interaction may challenge the scalability of previous 
schemes. We further propose a quantum computation scheme. In this scheme, by using appropriate 
encoding method, we effectively negate influence of the next-nearest-neighbor interaction in order 
to improve the precision of quantum gates. We also discuss the feasibility of this scheme in 1- 
D Josephson charge qubit array system. This work may offer improvement in scalable quantum 
computing. 
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I. INTRODUCTION 



One of the critical problems in realizing scalable quantum computation (QC) is perform- 
ing two qubits gates, which implies that the couplings between qubits are variable functions 
subject to external control. In many physical systems, this requirement is not easy to be 
satisfied. Recently various "always on" QC schemes have been proposed to solve this prob- 



lem. Zhou et al. suggest encoding logical qubits in interaction free subspace (IFS) 
while Benjamin et al. suggest changing interaction type between qubits from non-diagonal 
Heisenberg type to diagonal Ising type by tuning the Zeeman energy splits of individual 
qubits Pi . Schemes implementing above ideas into AFM spin ring and optical systems have 
already been developed 

These above elegant proposals mainly concern a 1-D spin-1/2 chain with couplings be- 
tween neighboring qubits. This model is a good analogue to many candidate scalable QC 
implementations. But in realistic systems such as quantum dot and optical lattice, there is 
not only the nearest-neighbor interaction, but also the next-nearest-neighbor, or even longer 
range interactions existing. In this paper we study effect of the non-nearest-neighbor cou- 
pling in 1-D spin-1/2 chain model. Since the strength of the long-range coupling is often 
much smaller than that of the nearest-neighbor coupling, in many previous schemes its ef- 
fect was ignored. But this omission results deviation in performing quantum gates. In this 
paper first we estimate the deviation of quantum gates induced by the omitted long-range 
interaction. Results imply that this deviation depends on not only the long-range interac- 
tion strength but also the length of the spin chain. This means that though the long-range 
interaction strength is very small, the induced deviation can "cumulate" on the whole 1-D 
spin-1/2 chain, hence scalability of previous QC scheme is restricted. Following this estima- 
tion, we consider how to suppress this unwanted effect. We propose a QC scheme in which 
by using proper encoding methods, influence of the perpetual next-nearest-neighbor inter- 
action is effectively negated, hence the precision of quantum gates updated. Compared with 
previous schemes in which the long-range coupling was omitted, our scheme does not cause 
the speed of quantum gates slow evidently. We also discuss the feasibility of this scheme in 
1-D Josephson charge qubit array system. 

This paper is organized as follow: In the second section we consider a 1-D spin-1/2 chain 
model. In the recent "always-on" QC schemes , long-range interaction terms are always 
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neglected. Here, we estimate the quantum gates deviation induced by an omitted perpetual 
next-nearest-neighbor Ising type interaction term. In the third section we consider how to 
suppress influence of the untunable long-range coupling. We use "blockade spin" encoding 
method to effectively neutralize unwanted influence of the perpetual next-nearest-neighbor 
interaction. Based on this encoding method we propose a QC scheme. In the fourth section 
we study the feasibility of this scheme in 1-D Josephson junction charge qubit array system. 
Several potential generalizations are suggested before conclusion. 



II. INFLUENCE OF THE PERMANENT LONG-RANGE INTERACTION 

We start with a 1-D spin-1/2 chain consisting of 2n + 1 spins, with tunable XX Z inter- 

n 

action between neighboring spin p. The Hamiltonian is: 



Hi deal = Hs + Hj (1) 



where 



2n+l 2n+l 



^= E^^= (2) 

i=l i=l 
2n 2n 

i=l i=l 

We assume that only values of 5*, B\, and Ji,i+i are tunable, while Ji remains constant. 
Various systems including single electron arrays, optical lattice, quantum dot, and Josephson 
junction array could be described by this "popular" Hamiltonian, while methods developed 
for the XX Z type exchange interaction can be easily generalized to other types of exchange 
interaction including XY type and Heisenberg type interaction j^. 

But in realistic systems there is not only the nearest-neighbor coupling, but also residual 
non-nearest-neighbor couplings existing. The interaction terms between non-neighboring 
qubits may origin from residual wave functions overlap or long-range Coulomb interaction. 
Without loss of generality, we may flrst estimate the influence from the next-nearest-neighbor 
coupling. We could assume that a permanent next-nearest-neighbor interaction term is 
ignored in Eq. ^ 

2n-l 

Hl = J2Y. (4) 



i=l 
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Thus though the theoretical evolution of the system is governed by Hjdeai, the realistic 
evolution is governed by the realistic Hamiltonian Hfi = H ideal + Hl. When we perform 



quantum gates 



concerned QlP ^ 



bllowing previous schemes in which only the nearest-neighbor interaction is 
deviation of the realistic evolution from the ideal expectation is induced 



f. H 



Let us follow schemes developed in ref. |3| as an example. For untunable ^ Jiafaf^^ term 

1=1 

in Hi, methods of freezing "blockade" spins in definite states between logical qubits have 
been proposed to negate the influence of this always-on coupling on single qubit operations. 
As shown in Fig. 1: In the spin- 1/2 chain, only the even spins (the hexagonal ones) are 
chosen as qubits while the odd spins (the rounded ones) are used as "blockade". When 
performing single qubit operations, we would make the blockade spins "frozen" in states 

2n 

|0) and |1) as shown in Fig. 1, thus the influence of the perpetual X] Ji^t^t+i terms on 

1=1 

the even spins is effectively neutralized because ^iO'2j(o"2i_i + ^^li+i) — ^^r any i. So 
there are n qubits and n + 1 "blockade spins" in the whole chain, where the ith qubit is 
encoded by the 2ith spin. The single qubit gates on the zth qubit are realized by tuning 
the effective magnetic field on the 2zth spin, while two qubits gates between the ith qubit 
and the {i + l)th qubit could be established by tuning the inter-spin exchange interaction 
strength J2i,2i+i and J2i+i,2i+2- 

Spin 1 23456 789 10 

2n 

FIG. 1: For the always-on ^ Jiu^afj^i type inter-spin coupling, qubits will suffer continuous 

i=l 

phase gates with their neighbors. We may only use the even spins (the hexagonal ones) to encode 
information while use the odd spins (the rounded ones) as blockade spin. When performing single 
qubit operation, the blockade spins are "frozen" in definite states |0) or |1) in order to negate 
influence of the perpetual Ising interaction. 

When we perform quantum gates, the theoretical expected unitary evolution of the whole 
chain is U = exp{—itHideai) while the realistic evolution is V = exp[—it{Hideai + Hl)]. 
We could estimate the deviation \\U — V\\ in single and multi qubit operations, where for 
convenience we use the definition of spectral norm, that is, the norm of an operator O is 
defined as the square root of the maximum eigenvalue of 0^0. As shown in appendix. 
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A, when we perform a certain quantum gate, the deviation induced by the omitted next- 
nearest-neighbor interaction is at least of the order \exp{—iJ2tn} — 1|, t being the required 
time for performing this gate. For t very small, we could estimate the deviation speed: 

j^\\U-V\\^0{nJ2) (5) 

This result implies that, though J2 is very small, deviation induced by Hi could become 
very large because this deviation depends on the length of the chain. This deviation could 
increase with the chain become longer. Thus omitting the long-range interaction effect may 
restrict the scalability of previous schemes. Below we consider how to handle this problem 
by using proper encoding methods. 



III. USING ENCODING SCHEMES TO SUPPRESS THE INFLUENCE OF 
LONG-RANGE INTERACTION 

Let us illustrate our idea intuitively. Following the previous section we start with a 1-D 
spin-1/2 chain with tunable nearest-neighbor XX Z interaction and next-nearest-neighbor 
Ising interaction. The Hamiltonian reads: 

Hm = Hs + Hj + Hl (6) 

Where Hs, Hj, and Hi are described by Eqs. 121111 Here, we further set Bl = for any i in 
the whole quantum information process, i. e. there is no cr^ terms in Hs- 

As sketched in Fig. 2, we encode one logical qubit by two physical spins, using the spin 
states |01) and |10) as logical |0) and Thus the tunable XY interaction between two 
spins plays the role of rotation in the 2-D logical Hilbert space 0|. 

Similar to the "blockade spin" methods used in ref. 0], our intuitive idea is freezing 
two "blockade spins" between each two logical qubits to negate the influence of the nearest- 
neighbor and next-nearest-neighbor always-on Ising type interaction. As shown in Fig. 2: 
The hexagonal ones are spins used to encode information while the octagonal ones are used 
as blockade. We use the spins 3 and 4 as one logical qubit while spins 7 and 8 as another. 
When performing single qubit operations, we would make the blockade spins (spins 1, 2, 5, 6, 
9, and 10) all "frozen" in state |0), thus the influence of permanent Ising type interaction on 
single logical qubit can be effectively negated: Since spin 3-4 are in Hilbert subspace spanned 
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by spin states |01) and |10), we simply calculate the influence of Ising type interactions on 
the logical qubit and get: 



(7) 



Spin 1 23456 789 10 




FIG. 2: We use two spins to encode one logic qubit, that is, the two spins 3-4 as one qubit and 
spins 7-8 as another. Two "blockade spins" both in |0) state are placed between two logical qubits 
in order to negate the influence of next-nearest-neighbor interaction. 

Now we show how to perform universal quantum gates, that is, the single qubit 
rotations, the single qubit rotations, and CPHASE gates between two qubits. As shown 
in Fig. 2, performing the single qubits rotations on qubit encoded by spin 3 and 4 is 
easy by tuning 73^4. We further note a trivial fact mentioned by ref. [l2| that the single 
qubit o"^ rotation can be constructed by the CPHASE gate and single qubit a^' rotation. 
So the central problem becomes performing CPHASE gate between two logical qubits. Our 
main idea to achieve this goal is, by tuning the exchange interaction J4 5 and Jg 7, we could 
perform a CPHASE gate between logical qubits encoded by spin 3-4 and spin 7-8, adding a 
phase on one of the four logical qubits states while remain the other three states unchanged. 

We separate the spins 3, 4 and 5 as one party and spins 6, 7 and 8 as the other. We 
set -B* = for any i in the whole process of performing two qubit gate. Starting with the 
four possible initial states {|100) , |010)}3 4 5 {|010) , |001)}g 7 g, initially we tune aU the 
strength of nearest-neighbor XY interaction to zero: Ji,j+i = for any i. Thus the the four 
possible logical qubits states have equal static energy and we deflne this static energy as 
energy zero point. As mentioned above, in performing two qubit gate, what we may tune is 
just J4,5 and Je,?, so we can reduce the Hamiltonian Hm in Eq. El into a 9-D Hilbert space 



The whole Hamiltonian Hm becomes a function of tunable J4_5 and Jej: = 
Hm{J4,5, Jgj)- Below we use the label \abcdef) to label the quantum states of the six 
spins from spin 3 to spin 8, the flrst a for spin 3, the second b for spin 4 the last / for 



S = span{{\m) , 1010) , |001)}3_4_5 (g) {|001) , |010) , \m},j^s} 



(8) 
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spin 8. For example, 1 100010) labels the state in which the spin 3 and spin 7 are in state |1) 
while the spin 4, 5, 6 and 8 are in state |0). 

Indeed, S in Eq. IHlcan be reduced to four Hilbert subspaces: The first one is a 1-D trivial 
subspace 5*1 spanned by single state 1 100001). The second is 2-D 



^2 = span{{ 1 100)} 



3,4,5 



{|oio),|ioo)W 



(9) 



Under basis {|100010) , |100100)} Hamiltonian Hm{J4,5, Jqj) can be reduced to: 



-f^2(<^4,55 Jqj) - 

The third is similar to the second: 

^3 = span{{|010),|001)} 



2J6,7 
2J6,7 



(10) 



3,4,5 



{|001)}6.r,8} 



(11) 



Under basis {|010001) , |001001)} the reduced Hamiltonian is 



2J4,5 

2J45 



(12) 



The fourth is 4-D 



^4 = {|010) , |001)}3^,5 (g) {|010) , |100)}g^7g 
Under basis {|010010) , |010100) , |001010) , |001100)} the reduced Hamiltonian is 



(13) 



i?4(<^4,55 Ji 



6,7; 



2J6,7 2J4,5 

2J6,7 -4J2 2J4,5 

2J4,5 -4J2 2J6,7 

2J45 2J6 7 -4Ji 



(14) 



The form of H^lJ^^^, Jqj) is quite similar to the NMR type Hamiltonian: If we set 
spa?7.{|010) , |001)}3 4 5 as the left "qubit" and span {\010) , |100)}g7g as the right, we could 
see that 74^5 and Jqj play the role of tunable local X operation for individual 2-D Hilbert 
space, and the perpetual Ising type interaction induce an untunable term of the type 
a(cr£ + cr|j.) + &cr£cr|j. in the Hamiltonian. Unfortunately, techniques developed for NMR 
such as refocusing usually can not be employed in other systems (especially many solid state 
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systems including quantum dot and superconducting circuits) because the assumption of 
fast, strong pulse can not be valid, so we choose an alternative way to perform CPHASE 
gate. 

Starting with the four possible initial states with degenerate static energy, 
{|100) , |010)}3 ^ (g) {|010) , |001)}g in the first step wc tune only 74^5 while set Jgj = 0. 
The XXZ interaction between spin 4 and 5 keeps states 1 100001) and 1 100010) unchanged. 

In space Ss taking Bloch sphere representation we see the transformation induced by J4,5 
is a rotation about X axis. But in space S4 the induced rotation is more complex. We 
note that since our initial state in 5*4 is only 1 010010), the induced transformation in S4, is 
finally reduced to a 2-D subspace S'^ = span{{\010) , |001)}3 ^ 5 {|010)}g g}. Under basis 
{ 1 010010) , 1 001010)} the reduced Hamiltonian is 



^4(^4,5,0) 



2J4,5 
2J4,5 -4J2 



(15) 



The induced transformation is a rotation about an axis on X — Z plane because the static 
energy of jOlOOlO) and jOOlOlO) are shghtly different due to the long-range interaction Hl. 
So in 5*4, with tunable 74,5 acting as u^-operation, we can perform a transformation 

exp {—i7ia^/2}, which is exactly a tt rotation around the X-axis in Bloch sphere represen- 
tation. We assume that the maximum value of 74^5 we could tune to is Xi/2. We set 



parameter 9 and X2 as cos 9 — 



= , sin^ = 



2J2 



sm29 



2J2 



Then we define rotation (we set h— 1) 



= , cos 2^ = 



X2 



y/xlM^'- 



f, and 



R{X) = exp 



iH^{X,0)tt/2JX^ + {2J2 



Thus we have 



R{Xi) = i 



sin 9 cos 9 
cos 9 — sin ^ 



R{X2) = i 

And the combined rotation in subspace 5*4 is what we want: 



sin 2^ cos 2^ 
cos 2^ -sin 2^ 



(16) 



(17) 



(18) 



R{Xi)R{X2)R{Xi) =exp{-i7r(7^/2} 



(19) 



8 



These above tuning of J4^5 at last form a unitary transformation Ui in the 9-D space 
S: It is nontrivial only in space ^3 and ^4: 1 100001) and 1 100010) remain unchanged; 
|010010) is transformed exactly into |001010); |010001) is transformed into superposition 
of |010001) and |001001). As shown in Eqs. ITBM191 the required time for this step is 
Ti = n/^Xf + {2J2f + n/2^l^cf+~{2J2f. We could see that Ti mainly depends the 
maximal XY interaction strength we could have. We further note an important fact that if 
we re-perform the tuning of 74^5 with inverse strength, i.e. Xi to —Xi, X2 to —X2, we can 
get the inverse operation of f/i in S. 

In the second step we tune 74,5 to zero but come to control Jqj. Quite similar to the 
previous step, this exchange interaction keeps states in Si and 83 unchanged while induce 
transformations in 5*2 and 5*4. In 5*2 under Bloch sphere representation the induced trans- 
formation is a rotation about X axis, while in subspace S4, since the initial state |010010) 
is exactly transformed into state |001010) by the first step, the induced transformation by 
Jqj in 5*4 is restricted to a 2-D subspace 5*4 = {|001)}3 4 5 {|010) , |100)}Q7g. Under basis 
(lOOlOlO) , lOOllOO)} the reduced Hamiltonian is 

The form of H'^{0, Jqj) is quite similar to that of the previous 1/4(74 5, 0), so we can perform 
another unitary transformation U2 just similar to the first step which implement a tt rotation 
around the X-axis in space 5*4, transforming state |001010) to |001100). 

After the above two steps, we review the intermediate states we get: 1 100001) remain 
unchanged; 1 100010) is changed into the superposition of 1 100010) and 1 100100); |010001) is 
changed into the superposition of |010001) and |001001); |010010) is changed into |001100). 
The previous three intermediate states are degenerate under Hamiltonian Hm (0, 0) but the 
last state has nonzero static energy — 4Ji. Therefore, in the third step we tune off all 
exchange coupling for a time interval r. In this period the state |001100) experience an 
additional phase cp due to its non-zero static energy. In the last step we can perform the 
inverse operation of U2 and Ui to transform the four intermediate states back to the initial 
four states. 

After all the above operations we have performed a CPHASE gate U {(f) between two 
qubits, adding a controllable phase (p on the state |010010) while remaining other states 



-4J2 2J6,7 
2Ja7 -4Ji 



(20) 
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unchanged. As we mentioned before, we use spin states |01) for spin 3-4 and spin 7-8 as 
logical |0) and |10) as logical |1), in this 4-D representation the gate U (y?) we obtain is 



10 
e^"^ 
10 
1 

After implementation of the CPHASE gate U, we note that 



(21) 
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thus with the single qubit cr^ rotation of the second qubit and CPHASE gate between the 
first and the second qubit, the single qubit rotation of the first qubit is obtained. 

In summary of this section, we have demonstrated all the required universal gates for 
QC. In this scheme, by using appropriate encoding method, influence of the next-nearest- 
neighbor interaction is effectively ruled out, thus the precision of quantum gates updated. 
Besides, the speed of the CPHASE gate mainly depends on the strength of nearest-neighbor 
exchange interaction, it is not restricted by the small value of J2. The quantum information 
speed of this scheme is in the same level with that of previous schemes in which the next- 
nearest-neighbor interaction was neglected. 



IV. A POTENTIAL PHYSICAL REALIZATION: JOSEPHSON JUNCTION 
CHARGE QUBIT 

Now we consider the application of our scheme to realistic systems. We consider the long- 
range interaction in Josephson charge qubit system as an example. The typical Josephson- 
junction charge qubit is shown in Fig. 3 1^: It consists of a small superconducting island with 
n excess Cooper pairs, connected by a tunnel junction with capacitance Cj and Josephson 
coupling energy Ej to a superconducting electrode. A control gate voltage Vg is coupled to 
the system via the gate capacitor Cg. The Hamiltonian of the Cooper pairs box (CPB) is 

H Josephson = Ec{n - Ugf + Ej COS $ (22) 

where E^ = (2e)^/2 (C^, + Cj) is the charge energy, Ug = CgVg/2 is the gate charge, and 
$ is the conjugate variable to n. When E^ ^ Ej, by choosing ng close to the degeneracy 
point Ug = 1/2, only the states with and 1 Cooper pairs play a role while all other charge 
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states in much higher energy level can be ignored. In this case the CPB can be reduced 
to an effective two-state quantum system. A further step is replacing the single Josephson 
junction by two identical junctions in a loop configuration in order to gain tunable tunnelling 



amplitude 



By making replacement n = {1 + a^) /2, the effective Hamiltonian can be 



written in the spin-1/2 notation as 

H Josephson 



(23) 



where the state with Cooper pairs corresponds to the spin state ||) and 1 Cooper pairs 
to It). and B^ are the effective magnetic fields which are controlled by the biased gate 
voltage and frustrated magnetic flux. 



E C 









(§) 

















FIG. 3: A Josephson charge qubit formed by a superconducting single-charge box and a tunnel 
junction. The single Josephson junction is replaced by a flux-frustrated SQUID so that the effective 
Josephson coupling is tunable. 

For coupling two CPBs, the direct capacitance coupling is most intrinsic, but its 

drawback is also obvious, that is, the coupling strength induced by connective capacitances 
is untunable. For a system consist of M CPBs coupled with each other by capacitances, the 
static charge energy can be written as |Q 



H, 



c 



{2ef 



nC rt 



(24) 



where rf = {rii — n^i, ^2 — ^^2, ''^3 — ''^gs, — ngM) is the charge number vector of the 

M CPBs, and C is the capacitance matrix of the system whose diagonal term Ci^i equals to 
the sum of capacitance around the ith CPB, and non diagonal term Cij corresponds to the 
connective capacitance between CPB i and j. 

A schematic plot of an array of N capacitively coupled CPBs is shown in Fig. 4. The 
CPBs have Josephson energies Eji and capacitances Cji. Each CPB is connected to the 
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control gate voltages Vgi via a gate capacitance Cgi. The ith intermediate CPB is connected 
to its neighboring (z ± l)th CPBs via the connective capacitors Cci,i±i- We assume that all 



CPBs are identical with the same capacitances: C, 



9« 



Cg, Cji 



Cj, and all the coupling 



capacitance Cd^i+i are equal to Cc- So we have a tridiagonal capacitance matrix for this 
system. For the intermediate qubits of the array, 



Ci,j = Co[Si,j{l + 2e) - 5i,j±ie], l<i,j <N 



(25) 



where Cq = Cg + Cj and e = Cc/Cq. For qubits on the edge of the array a small correction 
is needed: Ci^i = Cn,n = Cq{\ + e) while Ci_2 = Cn^i^n = —Cot. 






F C 



F C 



FIG. 4: A schematic illustration of array of capacitively coupled CPBs. 

Since C is a tridiagonal matrix, has nonzero matrix elements on the second, third 
and other diagonals which characterize the capacitance induced Coulomb interaction be- 
tween different CPBs. If e ^ 1, the off-diagonal elements of decay exponentially as 



C ■ ■ ~ 



. So the influence of the long-range interaction can be reduced by taking 



the coupling capacitances Cc much smaller than the on-site capacitances Cq. Again by mak- 
ing replacement Ui = {1 + al) /2, we get that the interaction term C^jijii — ngi){nj — Ugj) 
provides always-on Ising type interaction between CPB i and j. 

If we go further, replacing the coupling capacitance by SQUID, we can have tunable XY 
type exchange interaction between neighboring qubit besides perpetual Ising type interaction 
(Fig. 5) 13|14|l3- 

Due to flux quantization the phase across the coupling SQUID is 
$1 — $2 + ct; being a constant controlled by the frustrated flux. If we tune a to zero, the 
coupling term cos($i — $2) thus induce XY type exchange interaction crJ'cTg + CTiCrf- 

Thus we can see the correspondence between the theoretical Hamiltonian in Eq. IHl and 
the reahstic physical system: the tunable SQUIDs of single qubits induce tunable terms 
in the single qubit part of Hamiltonian; all the bias voltages are biased on the degeneracy 
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V. 





V. 



«2 



m 



FIG. 5: A schematic plot of two CPBs coupled by a SQUID. 

point Vgi = 1/2 in order to prevent the qubits from the 1/f noise effect implying 
that there is no term in the single qubit part of Hamiltonian; the exponentially decay 
capacitance coupling corresponds to the nearest Ising interaction part in Hj and the next- 
nearest-neighbor Ising interaction part Hl, while the tunable SQUID coupling corresponds 
to the tunable XY interaction part in Hj. We also see that the required performances in 
the above scheme just correspond to tuning external magnetic field frustrated in the SQUID 
loops. 

By taking e <^ 1 the effect of the long-range interaction in Josephson charge qubit array 
system can be reduced but can never be negated. Moreover, the speed of two qubits operation 
depends on e, decreasing the value of e would slow the speed of two qubits gate. Besides, 
making e smaller and smaller may be highly challenging in experimental realization. So we 
may choose an alternative way to handle the problem of always-on non-nearest-neighbor 
coupling. The main idea is that, based on the decay property of the long-range interaction, 
we can take some lower order terms of long-range interaction into account while taking the 
higher order terms as random noise. Thus we may use the encoding schemes discussed in 
Sec. Ill to negate the influence of the few lowest order couplings. 

Schemes of using capacitively coupled Josephson array to perform QC have been proposed 
Our scheme offer an alternative idea to handle the long-range interaction problem. 
The only parameter required to tune is the flux frustrated in SQUID loops, while the gate 
voltages are frozen on the degeneracy point, which prevents qubits from severe dephasing. 

It should be noted that, in our scheme although the next-nearest-neighbor interaction 
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effect is entirely negated, tfie fiigfier order long-range couplings still work. A natural gen- 
eralization is using 3 spins as one "blockade" to negate the influence of the third order 
long-range interaction (first order being the nearest-neighbor coupling while second order 
being the next-nearest-neighbor coupling), or even m spins as one "blockade" to negate the 
influence of the mth order long-range interaction. Following analysis similar to the estima- 
tion in appendix. A, we could conclude that, if we have negated the influence of the mth 
order interaction, the speed of deviation from ideal unitary evolution will be 0{nJm+i), Jm+i 
being the characteristic interaction strength of the (m-|- l)th order couphng. In 1-D Joseph- 
son charge qubit array system the strength of the long-range interaction have exponentially 
decay property, this means that, for a certain n, we could always make nJ^+i small enough. 
So we could apply the generalization of our QC scheme to this physical system to suppress 
the speed of deviation induced by long-range interaction into some tolerable domain. An- 
other potential extension is generalizing the nearest-neighbor coupling in our scheme from 
XX Z type to various other types including XY type and Heisenberg type. We could also 
translate our idea of handing long-range interaction in this paper to other implementations. 



V. CONCLUSION 



In conclusion, in this paper we have studied the non-nearest-neighbor interaction effect in 
1-D spin-1/2 chain model. We prove that the quantum gate deviation induced by the long- 
range interaction may challenge the scalability of quantum computing. We further propose a 
QC scheme in order to suppress influence of long-range interaction. In this scheme, using ap- 
propriate encoding method, we effectively neutralize influence of the next-nearest-neighbor 
interaction, thus the precision of quantum gates updated. The quantum information speed of 
this scheme is in the same level with that of previous schemes in which the long-range inter- 
action strength was ignored. We also discuss the feasibility of the scheme in 1-D Josephson 
charge qubit array system. This scheme may offer improvement in dealing with systematic 
errors in scalable quantum computing. 
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APPENDIX A: ESTIMATION OF DEVIATION INDUCED BY LONG-RANGE 
INTERACTION 

Let us follow schemes developed in ref. j^. As shown in Fig. 1: In the spin- 1/2 chain, 
only the even spins (the hexagonal ones) are chosen as qubits while the odd spins (the 
rounded ones) are used as "blockade" . 

When we perform quantum gates following schemes developed in ref. the theo- 

retical expected unitary evolution is U = exp{—itHjdeai) while the realistic evolution is 
V = exp[—it{Hjdeai + Hji)], We could estimate the deviation ||f/ — y|| in single and multi 
qubit operations: 

(1) Idle. When the whole chain is in idle status, the single qubit part Hs and the exchange 
interaction term in Hj are all tuned off, the n + 1 odd spins used as "blockade spins" are 
"frozen" in definite states |0) or |1). Since the permanent nearest-neighbor and next-nearest- 
neighbor Ising type interaction do not transfer energy, we could reduce our discussion in a 
2"-D Hilbert space H which is direct product of all the n qubits. The expected evolution is 
Uidie = I- But the realistic evolution is 

Vldle = exp{-iJ2t[-n + fT2i^^2i+2]} (Al) 

i=l 

(Tj, (t|, and a J being Pauli operators of the jth spin. Eigenvalues of Vjdie varys from exp{iJ2t} 
to exp{—iJ2t[—2n + 1]}. exp{iJ2t} corresponds to states that the n qubits are all in state 
|0) or all in state |1), while exp{—iJ2t[—2n + 1]} corresponds to states that for any i the 
2zth spin and the 2{i + l)th spin are in opposite states. We could choose a proper energy 
zero point for realistic evolution process, this means we could add a proper phase factor on 
Vidie '■ Vidie G^'^Vidie- But uo matter how we choose the energy zero point, we could prove 
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that, for any </?, either 

\e''^ exp{-iJ2t[-2n + 1]} - l| > \exp{-iJ2t[n - 1]} - 1| (A2) 

or 

\e''^ exp{iJ2t} - l| > \ex.p{-iJ2t[n - 1]} - 1| (A3) 
is vahd. So we could give an estimation of deviation for idle status: 

\\U - VI,,, > \eM-iMn - 1]} - 1| (A4) 

(2) Performing rotations. When performing rotations on the 2ith spin, parts for 
other spins in Hs and all Ji,i+i for any i are tuned off, the n + 1 odd spins used as "blockade 
spins" are "frozen" in definite states |0) or |1) as shown in Fig. 1. We could estimate the 

deviation \\Uz — Vz\\ in a 2"^^-D Hilbert space Hz- Hz is defined as below: Hz is a subspace 
of H. It is the direct product of all the n qubits except the ?'th, and for any state in Hz, the 
ith qubit is in state |0). The expected evolution in H^ is Uz — exp[—itBf]. But the realistic 
evolution is 

n-l 

Vz^eM-itiBl' -nJ2 + J2a'2i^2 + J2<+2 + J2 J2 ^2k^2k+2]} (A5) 

k=l,k^i—l,k^i 

Thus the estimation of deviation in H^ is quite similar to previous idle status: 

\\U - V\\f > \exp{-iJ2t[n - 1]} - 1| 

We note that the norm of operator U — V in the whole space H could not be smaller than 
the norm of operator U — V reduced in a subspace H^ oi H. So we give an estimation of 
deviation for performing gate: 

\\U - V\\^ > \exp{-iJ2t[n - 1]} - 1| (A6) 

(3) Performing rotations. Discussion of performing rotations on the ith qubit is 
similar to previous discussion of performing rotations on the ith qubit. Similarly we 
could define a subspace H^ of H: It is the direct product of qubits except the {i — l)th, 
ith, and (i + l)th qubits, and for any state in H^, the (i — l)th qubit is in state |0) , the 
ith qubit is in state (|0) + |l))/-\/2, and the {i + l)th qubit is in state |1) . So we have 
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^2C"|j_20"2i + ^20"2j'^2j+2 = 0- The cxpected evolution in this subspace is Ux = exp[—itBl^], 
while the realistic evolution is 

i-2 n-1 

K = exp{-it[Bl' - nJ2 + J2(yl(i^2) - ^2C^2(i+2) + + (A7) 

k=l k=i+l 

Eq. IA7I is similar to Eq. IA5I We give the estimation: 

\\U - V\\.^ > \exp{-iJ2t[n - 3]} - 1| (A8) 

(4) Performing inter-quhit gate. Inter-qubit gate between the ith qubit and the {i + l)th 
qubit is achieved by tuning the inter-spin exchange interaction strength J2i,2i+i and J2j+i,2i+2 
while other exchange interaction terms are all tuned off. Suppose in idle status the blockade 
spin 2i is in state |0), we could estimate the deviation in a 2""^ Hilbert space Hint- Hint is 
defined as follow: Hint is a subspace of H; It is the direct product of all qubits except the ith. 
and {i + l)th qubits. For any state in Hint-, the ii\i and the {i + l)th qubits are both in state 
|0). Obviously the ideal expected evolution in this space is /, but the realistic evolution is 

i-2 n-1 

Vint = exp{-it[-{n - 1) J2 - ^2^2(^-1) - ^2^2(^+2) + ^ Yl )^2fe<+2]} (A9) 

k=l k=i+2 

Quite similar to previous estimations we give: 

\\U - Vl^, > \exp{-tJ2t[n - 2]} - 1| (AlO) 
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